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LOW-DISCREPANCY SEQUENCES FOR PIECEWISE SMOOTH FUNCTIONS 
ON THE TWO-DIMENSIONAL TORUS 


LUCA BRANDOLINI, LEONARDO COLZANI, GIACOMO GIGANTE, AND GIANCARLO TRAVAGLINI 


Abstract. We produce explicit low-discrepancy infinite sequences which can be used to 
approximate the integral of a smooth periodic function restricted to a convex domain with 
positive curvature in R^. The proof depends on simultaneous diophantine approximation 
and a general version of the Erdos-Turan inequality. 

Keywords: Koksma-Hlawka inequality, piecewise smooth functions, discrepancy, dio¬ 
phantine approximation, Erdos-Turan inequality. 


1. Introduction 


Let / be a suitable function on the tf-dimensional torus /l/, and let {t ( 

be a distribution of points on T'^. The quality of the approximation of Jjd f{t)dt by the 
Riemann sum T!j=i f (t U)) is a basic problem with applications in 2D or 3D com¬ 
puter graphics, and also with applications when d is large (and the curse of dimensionality 
appears). See e.g. a. Any bound of the form 



can be termed a Koksma-Hlawka type inequality, provided the RHS is a variation V (/) of 
the function / times a discrepancy D of the hnite set {t with respect 

to a reasonably simple family of subsets of T'^. 

The case = 1 is the amazingly simple Koksma inequality, where T is replaced by the 
unit interval, V (/) is the usual total variation and D ({t j is the *-discrepancy 


1 ^ 

77 E ^[O.a) (K;)) - « 


0<a<l ” j=i 


that is the discrepancy measured on the family of all intervals anchored at the origin. 

See |l3], ||9], HD, lfl5l . llT], i22\ . ll^ as general references. 

The term Koksma-Hlawka inequality properly refers to E. Hlawka’s generalization of 
Koksma inequality to several variables, where / is required to have bounded variation in 
the sense of Hardy and Krause. In one variable, many familiar bounded functions have 
bounded variation, but, in several variables, the Hardy-Krause condition cannot be applied 
to most functions with simple discontinuities. For example; the characteristic function 
of a polyhedron has bounded Hardy-Krause variation if and only if the polyhedron is a 
-dimensional interval. 

We recall some of the variants of the Koksma-Hlawka inequality which have appeared 
in the literature so far. F. Hickernell IMI has proposed Koksma-Hlawka type inequalities 
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for reproducing kernel Hilbert spaces. J. Dick ifTSll has used fractional calculus to prove 
a Koksma-Hlawka type inequality for functions with relaxed smoothness assumptions. G. 
Harman lfT9l has considered a geometric approach and measured the variation by counting 
the convex sets needed to describe super-level sets of the function /. In lb) the authors of 
the present paper proposed a Koksma-Hlawka type inequality especially tailored for sim- 
plices, while in fTl they have introduced a Koksma-Hlawka type inequality for piecewise 
smooth functions. Analogues of the above problem in more general settings can be found 
e.g. in m and 15]. 

K. Basu and A. Owen H] have recently produced low-discrepancy sequences for a tri¬ 
angle, where the discrepancy is the one considered in 0. In this paper we propose a 
sequence of points which gives low discrepancy in the sense of the Koksma-Hlawka type 
inequality in jT]. We first recall a particular two-dimensional case of the statement therein. 


Theorem 1 (Q). Let h[t) = f (t) Xa {t), where f is a smooth 7?-periodic function on . 
and Xq. A the characteristic function of a bounded Borel set in Let 

d^f 


V(/):=4||/||^,(t2)+2 
Let{t{j)}%,C. 


df 


dt\ 


+ 2 


df 


dtj 


L'(t2) 


H (t2) 

, for any s S (0,1)^ and for any x S R^ let 
I{s,x) =U^ez2([0>'5i] X [0,S2]-\-x-\-m) 


dt\ dti 


L‘(t2) 


and let 


( 1 ) D({t{j)}^^:= sup E (;)+'«) 


Then 


^E E h{t{j)+m)- f h{t) dt < V(/) D[{f(;)}^=i) • 


Observe that if a set K € R^ does not intersect any of its integer translates, then it can 
be thought of as a subset of T^, and in that case the expression 

^E E XKitij)+m)-\K\ 

compares the measure of K with the share of points in K of the collection obtained by 
projecting {f(y)}^i onto T^. It follows that the above theorem includes, but is slightly 
more general than the analogous theorem where not only the function / but also the set 
n and the point collection {f(y)} are in T^, and the quantity D({f(y)}^j) is just the 
discrepancy with respect to the intersection of Q. with all the rectangles in T^. 

We are therefore interested in choices of the set {t (y)}^i which give satisfactory upper 
bounds for the discrepancy ([TJ. 

An interesting result in this direction is due to J. Beck 0: for every positive integer N 
there is a collection of N points in the unit square with isotropic discrepancy (that is, the 
discrepancy with respect to all convex sets) bounded by cA^^/^log^A. Since the discrep¬ 
ancy O is smaller than the isotropic discrepancy. Beck’s result gives a sequence that can 
be used in the Koksma-Hlawka type inequality in Theorem [T] when Q. is convex. On the 
other hand. Beck’s construction is somewhat intricate, and is obtained partly by random 
and partly by deterministic methods. 
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A more explicit extensible construction comes from a result of H. Niederreiter (see 
||2^ or II 2 T] page 129 and page 132, Exercise 3.17]): if 1, a, j3 are algebraic linearly 
independent on Q, then the discrepancy of {{ja, with respect to all axis parallel 

rectangles contained in the unit square is bounded by . This immediately implies 

that the isotropic discrepancy of this sequence is bounded by (see ETl Theorem 

1.6, page 95]), an estimate that is far from Beck’s result. 

Our main result is the following. 

Theorem 2. Assume that a,li are real algebraic numbers and that 1,CC,I3 is a basis of a 
number field on Q of degree 3. For all integers j > 0, let t (j) = (ja,jp). Let LI be a 
convex domain contained in with boundary having strictly positive curvature. Then 
the discrepancy defined in 0 satisfies 

(2) <cA-2/3iogA. 

The above constant c depends on the minimum and the maximum of the curvature ofdLl, 
on its length, and on the numbers <X, fi. 

For example, one can take a = where ^ is a real root of a third degree 

irreducible polynomial in Z. 

In other words, Theorem|2]says that a regularity assumption on the convex set fi suffices 
for the sequence in Niederreiter’s result to improve Beck’s estimate log^ A. This can 
be obtained by estimating directly the discrepancy D{{t{j)}^^i), and avoiding the isotropic 
discrepancy. The main tool that will allow us to do it is a version of the Erdos-Turan 
inequality essentially contained in mi. 


2. Proofs and auxiliary results 


Let us begin by recalling the above mentioned general form of the Erdos-Turan inequal¬ 
ity. 


Theorem 3. There exists a positive function \j/ [u) on [0, -f 0 °) with rapid decay at infinity 
such that for every collection of points {f (y)}y=i C for every bounded Borel set D C 
and for every R>0, 

L XD{tU)+m)-\D\ 

j= 1 rnGlP- 


< 


Hr{0) 


E 

neJ? ,Q<\n\<R 


%{n)\ + HR{n) 


1 ^ 

-T 


^2mn-t{j) 


Here 

Hr{x) = \j/{Rdist{x,dD)), 
where dist is the Euclidean distance in R.'^. 


Proof Take a smooth radial function ot(^) supported in |^| < l/2andwith = 

1, and define 


A (x) = (1 + 1^ P) (,« , m) , 

K(y)dy'] [ K{y)dy 
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Since ^(^) = 0 if |^ | > 1, it follows from the Paley-Wiener theorem that K{x) is an entire 
function of exponential type smaller than 1 , positive with mean 1 , all its derivatives have 
rapid decay at infinity, and K{^) <1 for every ^ e If we set Kr (x) = R^K (Rx), then 
the functions 


A{x)= [ KR{y){xD{x-y)-HR{x-y))dy 

Jm< 

Six) = [ KR{y){xD{x-y)+HR{x-y))dy, 

Jm< 

are entire of exponential type smaller than R and 

^ (^) < Xd (^) ^ B{x), \B (x) — A (x)| < 4\j/ {Rdist{x,dD)/2) 
(see HD for the details). Periodization gives 

^ A (x + m) < ^ Xoix + m) < ^ B{x + m), 

and, by the Poisson summation formula, 

^ A(x + m)= ^ K{R-\)(xD{n)-HR{n)y^^‘--\ 

^ B{x + m)= ^ (2d (n)+//«(«)) 

m€’L^ 


are trigonometric polynomials of degree at most R. It now follows that 
E XD{t{j)+m)-\D\ 


- 4E E B{t{j)+m)-\D\ 

= i!iE E K{R-^n)(xD{n)+HR{n))e^-'^<^'>-\D\ 

j=lneZ‘l 

= HRiO)+ ^ K {R-\) (2d {n)+HR («)) ^ f 

n^Z^,0<\n\<R J=1 


< 


HRiO)\+ Y. (l 2 D(n)| + |//tt(n)|) 

,0<\n\<R 


1 N 

— y 


Similar estimates from below can be proved, if one uses A (x) instead of B (x). 


□ 


A second tool in the proof is the estimate of the Fourier transform of arcs of curves in 
The next two lemmas are well known (see e.g. lIZTl Chapter 8 ]). We recall the proof 
of the first one both to help the unfamiliar reader, and to emphasize its two-dimensional 
nature. 

In what follows, for any arc 7 we will denote with y{^) the Fourier transform of its 
arclength measure. 

Lemma 4. Let Q,be a convex set in with a boundary with non-vanishing curvature. 
Let ybe an arc ofdLl and tCmin >0 be the minimum of the curvature of y. Then for |^ | > 1, 
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the Fourier transform is bounded by 

Here £ is the length of the arc and c is a universal constant. 

Proof. Let r (t) be the parametrization of 7 with respect to arc length, so that 

7(^)= 

Jo 

For any ^ we have the trivial estimate 


Jo 


< 1 . 


Assume 7 ^ 0 and let 

^=pri 

where 1 1 ] | = 1 and p > 0. First consider the (at most) three intervals f, f and f where 
|r'(T)-T 7 | > By Van der Corput’s lemma, since |r'(T)- 77 | > 2^^/^ and the ex¬ 

pression r" (t) • 77 = — K'(t) V (t) • 77 changes sign at most once (here v (t) and K'(t) are 
respectively the outer normal and the curvature of 7 at a point r (t)), then 


i 




< 


Cl 


(i = 1,2,3). The constant ci is universal. If |r' (t) • 77 1 < 2 we have |v (t) • 77 1 > 2 
so that 

|r"(T)-77| = Jf(T)|v(T)-77| > K-min2^*/2. 

Thus, by Van der Corput’s lemma, for the at most three intervals Ji,J 2 and Jj, where 
V (t) • PI < 2 ^^/^, we have 




< 


Cl 


(K'minP) 


(7 = 1,2,3). Again, ci is a universal constant. Finally, 


/ 


-2mpr(z)-ri 


< min f 


3c I 


1/2 


3c2 


(K-minP) 


1/2 


When p > 1 , this gives 


f\-2Kipr(xyn^^ 

Jo 


/ 1 -I- JC 

<min f,c 


1 / 2 ' 


,1/2 


□ 


Lemma 5. The Fourier transform of the arclength measure on the segment 7 joining two 
points X and y in is 

In particular, calling 1= |x —y| and 9 = 'T'c have 

| 7 (<^)| <min (e. 


7t\^-e\ 
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Proof. This is just an explicit calculation. 


□ 


Before we proceed with the proof of Theorem |2] we need a few results on convex sets 
in K'^. Let us begin with some terminology. 


Definition 6. Let K be a non-empty compact convex subset (a “convex body”) o/R'^. The 
signed distance function 5 k is defined by 


~ _ J dist(x, dK) if X G K 

Ok{x) = I ^ ^ 

For any real number u, define 

K“ = |x G R.'^: 5k{x) > m| 

and 

Ki, = |x G : 5k{x) = m| 

The signed distance function is Lipschitz continuous with constant 1, and IVS/rl = 1 
almost everywhere (see iflTl Section 14.6]). 


Definition 7. Let B be the closed unit ball centered at the origin. IfK is a convex body in 
R"^, then the outer parallel body of K at distance r is defined as the Minkowski sum of K 
and rB, 

K-\-rB = {x-\-y:x€K, |y| < r}, 

while the inner parallel body of K at distance r is defined as the Minkowski difference of 
K and rB, 

K ~ rB = {x '. X rB C. K}, 

Lemma 8. Let K be a convex body in R'^. 

(i) For any real number u, the set K'^ is the outer or the inner parallel body ofK at distance 
\u\, according to whether u is negative or positive, that is 

K" = KF\u\B,ifu<Q, 

K" = K^uB,ifu>Q. 

(ii) For any real number u, the set K" is convex (possibly empty). 

(Hi) IfM is a convex body too, then for every u > 0, 

{Mf\K)^ = (M„ r\K'*) U (M" r\Ku). 


Proof. Point (i) follows easily from the definitions, while the proof of (ii) can be found in 
||2^ Chapter 3]. As for point (iii), we sketch a proof, highlighting the main steps. First 
observe that (9(/T") =Ku and that (Mfl/T)" =M'‘r\K“ when m > 0. The thesis now follows 
after the observation that for any two compact sets A and B one has 

d{AnB) = {dAr\B)u(AndB). 


□ 


Lemma 9. Let K be a convex body in R'^ with boundary and let (Cmax be the maximum 
of all the principal curvatures ofdK. Finally, let 

r = r(A", tfmax) = {x '. (2K'niax) < 5k{x') < (2K'niax) }• 
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Then 5 k S Furthermore, the level set Ki, is whenever \u\ < * and its 

principal curvatures at a point x are given by 

l-uKj{y) 

where y is the unique point of dK such that dist(x,y) = \u\ and t^j{y) are the principal 
curvatures of dK at y. 

Proof This is essentially a reformulation of Lemmas 14.16 and 14.17 in El for the case 
of convex bodies. □ 


Let us now move back to the two-dimensional case. In the next two lemmas we estimate 
the Fourier transforms of the functions Xd and Hr in Theorem [2 for the specific type of 
sets D that one needs in the proof of Theorem |2] 


Lemma 10. Let FI be a convex body in with boundary with non-vanishing curvature 
and let KVnin cind tCmax be the minimum and the maximum of the curvature of dFI. Let I be 
a rectangle contained in a unit square with sides parallel to the axes, and call K = FI Cl 1. 
Then there exists a constant c depending only on KVnin such that for all R > and for 

every n = {n\,n 2 ) S with 0 < |n| < R, 

i'- / M 1 1 1 

\HR{n)\ < c--^+c 

\nfl^ 1 + |ni| 1 + \n2\ 

Here, Hr{x) is the function defined in Theorem\^by Hr[x) = y/(7?|5A-(x)|). Finally, there 
is a universal constant c > 0 such that for all R> 1, 


Proof By the coarea formula (see M Theorem 2, page 117]), since |V5 a:(x)| = 1 almost 
everywhere, 

HR{n)= [ \i/{R\5K{x)\)e^^^’’‘"dx 

= [ \j/{\Ru\) f e^^^‘^'"dxdu. 

i-oo Jk„ 

where Ku = {x : 5 k{ x) = u} as in the above Definition |6l and the integration on the level 
set K,, is with respect to the Hausdorff measure. Thus 


Hr (n) 

< / y/(R\u\)du sup 

/ e-^^‘^''dx 


J\u\<R-F2 |u|<ff-l/2 



+ 


I 

J\u 


\u\>R-'^P 


< — sup 

^ |«|<K-l/2 


WiR\u\) \Ku\du 


L 


e-^^‘^-'‘dx 


+ 


C2 


R 


10 ' 


The constant ci is just the integral of 21/ on [0, +oo), while C 2 depends on the rapid decay 
of 1 / and the slow growth of \Ku\ (recall that K" is convex and contained in a square of 
side 1 +2 |m|, and therefore the Hausdorff measure of Ku is smaller than 4(1 4-2 |m|)). In 
particular, ci and C 2 are universal constants and we immediately have that for any R>1 


Hr{0) 


< 


R' 


where c is a universal constant. 
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K 



Figure 1 . A convex body K — ^r\I and the relative set K^, with u <0. 
n has smooth boundary with non vanishing curvature and / is a rectan¬ 
gle. 

Now assume n ^ 0, < 1/ (2K'n,ax) and 0 <u Then, by the above Lemma 

0and Lemma|9] Ku consists of at most four smooth convex curves with curvature bounded 
below by tCmin, and at most four segments of length at most 1 parallel to the axes. By 
Lemma|4]and Lemma|5]this gives 



where the constant c depends only on the curvature KVnin- On the other hand, if < 

1 / (2K'niax) and if —< m < 0, then is composed by at most four smooth convex 
curves with curvature greater than or equal to 2ifniin/3, at most four segments parallel to 
the axes and of length at most 1, and at most eight arcs of circles of radius |m|. In order to 
better understand this, observe (see Figure[T]i that one can divide the complement of K into 
at most sixteen regions by taking the two normals to dK at each “vertex” of K (there are at 
most eight “vertices”). The part of Ku that intersects a region attached to a straight line is 
a parallel straight line of length at most 1. The part of Ku that intersects a region attached 
to a curve coming from (90 is a part of 0„. Finally, the part of Ku that intersects a region 
attached to a vertex of K is an arc of circle of radius \u\. 

It follows that 




































LOW-DISCREPANCY SEQUENCES 


9 


where the constant c depends only on the minimal curvature fCmin- Therefore, when 0 < 
\n\ < R we have 


Hr (n) 


11 ^11 1 
-"r |„|V2+^,5f;i + |„,|+Vo 

-^^^^i + lml i + kl' 


□ 


Lemma 11. Let Lite a convex body in with boundary with non-vanishing curvature 
and let fCmin be the minimum of the curvature of dLl. Let I be a rectangle contained in a 
unit square with sides parallel to the axes, and call K = Lid 1. Then there exists a constant 
c depending only on iCmin such that for every n = (ni,n2) S \ {(0,0)} 

, M 1 1 1 

\XKin)\ < c-^+c 

\n\^H 1 + |ni| 1 + |n 2 | 


Proof An application of the divergence theorem gives 


\XK{n)\ 


f e-^^‘''-^dx 


1 

JK 


JdK 2ni\n\ 


1 

2n\n\ 



Here v(;c) is the outer normal to dK at the point x. This oscillatory integral can be estimated 
by means of standard techniques. We include the details for the sake of completeness. The 
boundary of H n / is composed by at most four smooth convex curves with curvature 
bounded below by iCmin, coming from dPi, and at most four segments of length at most 1 
parallel to the axes, coming from dl. We therefore split the above integral into a sum of 
integrals over the components of dK described above. When integrating over a segment, 
the quantity v{x)-n/\n \ remains constant and an immediate application of Lemma|5]gives 
the estimate 

1 1 

*^1 + |ni| 1 + |n 2 | ’ 

with c a universal constant. Let us now estimate the integral over an arc of dPi, call it y. If 
r{'c) is a parametrization of y with respect to arclength, integration by parts gives 



Here K'(t) is the curvature of y at the point r(T) and Jq K(T)dT is the total curvature of 
y. Since y is an arc of dQ., the total curvature of y is smaller than the total curvature of 
do., that is In. The last inequality is just an immediate application of LemmalU where the 
constant c above depends only on the minimal curvature KVnin of dLl. 


































10 


L. BRANDOLINI, L. COLZANI, G. GIGANTE, AND G. TRAVAGLINI 


□ 


We are now ready to proceed with the proof of the main result of the paper. 

Proof of Theorem^ Let KVnin and tCmax be the minimum and the maximum of the curvature 
of do.. If we call mi,..., the lattice points for which the sets 

([0,ii] X [0,i2]+.r + m,') nil 

are nonempty, and let 

Ki = ([o,si] X [o,s2] +A+m,) nn, 

then of course 

^^€22 (([0>'5i] X [ 0 ,i 2 ]+x + m)nn) = vf.^^Ki. 

The number q is bounded by the maximum number of unit squares with integer ver¬ 
tices that intersect any given translate of H in This number is of course bounded 
by (diam (Q.) + 2)^ . We recall that we need a uniform estimate with respect to s and a. 



Figure 2. The intersection of a convex set D. with smooth boundary 
having non-vanishing curvature with the integer translates of a fixed rec¬ 
tangle. 


The sets Ki are as in Figure |2j at most four sides are parallel to the coordinate axes, 
while the curved parts come from dO.. The discrepancy 

1 ^ 

7=1 
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is clearly bounded by the sum of the discrepancies of the sets A', , 


L 


(=1 


L XKiitU)+m)-\Ki\ 


and we shall therefore study the discrepancy of a single piece Kj. Let us call K one such 
set. 

By the general form of the Erdos-Turan inequality in Theorem |3l the discrepancy of a 
single piece K is bounded by the quantity 


(3) 


HRiO) 


+ E (\XK{n)\ 

n^'L'^ ,Q<\n\<R 


HR{n) 


1 ^ 


We recall that/? > 0 is a number that we can choose at our convenience, Hr{x) = i^(/?|5a-(ji:)|) 
and V^(m) is a properly chosen function on [0, + 0 °) with rapid decay at infinity. 

The estimates of 2 a’(«) and //«(«) are contained in the above Lemmas fTOl and fTTl while 
the estimate of the exponential sums follows a standard argument. 


1 

N 


E 


;=i 




1 N 

— Y e2;Ei7« (a,j3) 


1 sin{nNn • (a,j3)) ^ 1 

N sin(;rn • (a,j3)) “ A^||n ■ (a,)3)|| ’ 

where ||m|| is the distance from u to the closest integer. 

Overall, the goal estimate Q becomes 


1 _ / 1 1 1 \ 1 
^ ^0<h<R V ^ 1 + |«l I 1 + H )N\\n- (a,^)|| ■ 


Observe now that 


log«log« , , 2'+'-12J+*-l 1 

■ ^ S S 2' 22 „,?2; II«1« + «2^ II 


+-Ei E 


i=0 2' „“2' II”'“I 


logff , 22+'-! 1 


7=0 2^' II«2/3||' 


Let US Study the sum I|22i0! + «2j3|| ^ first. By the celebrated result of W. 

M. Schmidt ll24l . see also ll25l Theorem 7C], since 1, oc, j3 are linearly independent on Q, 
for any £ > 0 there is a constant 7 > 0 such that for any n ^ 0, 


(4) 


||nia + n2/3|| > 


(l + |ni|)'+®(l + |n2|)‘+"' 

Then, following m, in any interval of the form 


[k-\)Y 


ky 


(1 +2''+')'+''(l +22+l)'+e ’ (1 +2''+')'+«(l +22+l)'+e 


where A: is a positive integer, there are at most two numbers of the form ||ni a + n2j31|, with 
2' <n\ < 2'+' and 2^ <n 2 < 2^+*. Indeed, assume by contradiction that there are three 
such numbers. Then for two of them, say ||ni a + 112 ^ || and ||mi a + 1712 ^ ||, the fractional 
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parts of «! a+ 112)3 and mia + m 2 ^ belong either to (0,1 /2] or to (1 jl, 1). Assume without 
loss of generality that they belong to (0,1 /2]. Then 

> |||nia+ 112)3 II - ||OTia + i«2)3||| 


(1+2<+1)1+£(i+22+1)1+'! 


= |iiia + ii2)3 — p — (i«ia + i«2)3 —q)\ 

> ||(i!i -mi)a + (ii2-'«2))3|| 

7 


> 


(l+2+i)i+c(l+22+1)1+'! ■ 

By the same type of argument, in the first interval 0, (i^ 2 '+i)i++i+ 22 +i)i+£ ) ’ 
points of the form ||iii a + nifi ||. It follows that 


there are no 


2 '+ 1 - 122+*-1 

L L ■ 

n\=2^ n2=2j 


Similarly, 


2- ll«l«l 


1 2+2 2{+;)(l+r) ^ 

I-^ <cV -^-< c2('+^)(i+®) (i + j). 

|iiia+ 112)3 II ky 


22+>-l 1 

II«2P|' 


n2=2J 


Finally, 


1 1 


0 <H</t(^ + l"ll) (l + l” 2 |) ll«-(«. 


logniogff , , iugj\ , iug/\ , 

<«E E 4+'+'+'’(i+2)+cE +'"+-'+cE 


\ogR 


logfi 


2=0 ;=0 
log« 


2 ' 22 


!=0 


J=0 


22 


<cY 2‘^R^logR + cR‘'logR < cR^^logR. 


=0 


Finally, we use the hypothesis that 1, a,)3 are a basis of a number field in Q. By a simple 
argument in number field theory, there is a constant rj such that for any n^Q, 

V 


||iiia+ 112)3 II > 


(max(|iii|,|ii 2 |)) 


2 ■ 


See for example 11251 Theorem 6F]. By a similar argument as before, this implies that 


Thus, 


2'+'-l 

L 

max{|ni|,|« 2 |)= 2 ' 


1 1 


1 


22' -2! 


(«,)3)ll 


< c ^ — < ci2^‘. 


/c=l 


logR 2'+*-l 1 1 

<"E E 


0<\n\<R\nf^^ ll«’(0!i^)ll !-=0 max(|«i|.|n 2 |)= 2 '' ll« ' («: 

log 2 ? 2 

< c £ :^i22'' < 

!=0 ^ ' 


Setting gives the desired estimate N log A, as long as A > 
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